High-speed Jet Formation after Solid Object Impact 
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A circular disc impacting on a water surface creates a remarkably vigorous jet. Upon impact an 
axisymmetric air cavity forms and eventually pinches off in a single point halfway down the cavity. 
Immediately after closure two fast sharp-pointed jets are observed shooting up- and downwards 
from the closure location, which by then has turned into a stagnation point surrounded by a locally 
hyperbolic flow pattern. This flow, however, is not the mechanism feeding the two jets. Using 
high-speed imaging and numerical simulations we show that jetting is fed by the local flow around 
the base of the jet, which is forced by the colliding cavity walls. Based on this insight, we then show 
how the analytical description of a collapsing void (using a line of sinks along the axis of symmetry) 
can be continued beyond the time of pinch-off to obtain a quantitative model for jet formation which 
is in good agreement with the numerical and experimental data. 

PACS numbers: 47.55.N-, 47.55.D-, 47.55. df, 47.11. Hj 



The most prominent phenomenon when a solid object 
hits a water surface is the high-speed jet shooting up- 
wards into the air. The basic sequence of events lead- 
ing to this jet has been studied since Worthington more 
than a century ago: After impact, the intruder creates 
an air-filled cavity in the liquid which due to hydrostatic 
pressure immediately starts to collapse, eventually lead- 
ing to the pinch-off of a large bubble. Two very thin 
jets are ejected up- and downwards from the pinch-off 
point which constitutes a finite-time singularity inten- 
sively studied in recent time P, [1, B S @1 ■ Such sin- 
gularities have been shown to lead to a hyperbolic flow 
pattern after collapse and thus to the formation of liquid 
jets 0,ii. 

As we show in the present work, however, the radial en- 
ergy focussing towards the singular pinch-off point alone 
is not sufficient to explain the extreme thinness of jets 
observed after the impact of a solid object, fnstead, this 
jet formation is shown to depend crucially on the kinetic 
energy contained in the entire collapsing wall of the cavity 
even far above the pinch-off singularity. 

This is in sharp contrast to jets observed in many other 
situations where narrow confining cavity walls are not 
present, e. g. for bubbles bursting on a free surface or 
near a solid wall [1, [13, [lH, impact of liquid droplets 
EEQEIll, wave focussing [13, [11 or jets induced 
by pressure waves [l^. fn all these cases jetting must 
thus be accomplished by a mechanism different from the 
one to be described in this paper. 

Our experimental setup consists of a circular disc with 
radius Rq that is pulled through a water surface with 
velocity Vb as described in [J. The velocity Vq is kept 
constant throughout the whole process. Global and local 
Reynolds and Weber numbers are fairly large as shown 
in [J] and therefore the only relevant control parameter is 
the Froude number, Fr= Vq / Rog with gravity g, which 
we choose to equal 5.1 (for Rq = 2cm and Vq = Im/s). 



Due to the large Reynolds numbers involved [23| and 
as during the highly transient process of jet formation 
only a negligible amount of vorticity is introduced into 
the system [lO[ we treat the problem as inviscid and irro- 
tational allowing for a potential flow description. Wc cm- 
ploy an axisymmetric boundary-integral technique which 
explicitly tracks the free surface. The topology change at 
pinch-off is implemented as follows: When the radial po- 
sition of the node closest to the axis becomes smaller than 
the local node distance, the two neighboring nodes are 
shifted to the axis, conserving their vertical position and 
their potential. Continuing the simulation, these nodes 
eventually form the tip of the top and bottom jets. These 
numerical simulations have shown excellent agreement 
with experiment for different impact geometries [1, [2^ 
and we verified carefully that our results are indepen- 
dent of numerical parameters such as node density and 
time stepping. All simulations include surface tension 
but neglect the influence of air. 

Figure [T] shows the pinch-off of the impact cavity and 
the subsequent formation of two thin jets. We use polar 
coordinates with z = at the pinch-off height and t = 
at the pinch-off moment. Velocity, length and time scales 
are normalized by Vq, -Ro, a-nd Tq = Rq/Vq, respectively. 

We set out to elucidate the precise mechanism which 
turns the horizontally collapsing cavity of Fig. [1] (a) into 
the thin vertical jets of Fig.[T](b) and (c). For this we fo- 
cus our attention on the dynamics of the upward jet base 
which we define as the local surface minimum illustrated 
in Fig. [21 ft is remarkable how the geometric confinement 
of the narrow cavity forces the jet to move upwards very 
fast while the widening of its base is restricted by the 
collapsing walls. We find that jet formation occurs on an 
extremely short time scale: the upwards jet grows above 
the initial quiescent surface in less than 1% of the total 
time after impact. 

These high speeds, however, are not due to a hyper- 
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FIG. 1: The free surface shape (black) and the disc position 
(red) from the simulation at pinch-off (a), at an intermediate 
time with the growing up- and downward jets (b) and at the 
instant when the downward jet hits the disc (c) . 



0.8 



0.6 



DC 




FIG. 2: The inset illustrates the position of the jet base (red 
diamond) at different instants of time. In contrast to other 
physical jetting situations the upward motion is much faster 
than the widening of the jet base. The main figure shows the 
upwards motion Zb{t) of the jet base predicted by the analyt- 
ical model (black line) which compares very favorably with 
the simulation (red crosses) and the experimental data (blue 
diamonds). (The slightly slower motion in the experiment 
can be attributed to an imperfect axisymmetry reducing the 
radial focussing effect.) The agreement between the model 
and the numerics is equally good for the widening of the base 
ri,{t). Here, due to optical difficulties no experimental data 
are available. 



bolic flow around the original pinch-off point as one could 
have expected based on suggested jetting mechanisms in 
other situations 0, B @ • Figure [3] demonstrates how the 
fluid here is not accelerated upwards continuously from 
the pinch-off singularity but instead acquires its large ver- 
tical momentum in a small zone located around the jet 
base: Since each horizontal cross-section of the axisym- 




FIG. 3: The vertical material acceleration = Dvz/Dt at 
t/To — 0.028 is confined to a small region around the jet base. 
The pinch-off location at (0,0) lies far too deep to influence the 
jetting process longer than in the first instants after pinch-off. 



metric cavity wall will keep flowing radially inwards even 
after pinch-off, eventually it must collide on the axis in 
a similar way as the original pinch-off. This creates an 
upward and downward acceleration, of which the upward 
acceleration feeds the jet. The downward (negative) ac- 
celeration below the jet base can clearly be observed in 
Fig. [21 It is thus essential to consider not only the sin- 
gularity itself but the continuous collapse of the entire 
cavity wall in any kind of theoretical modelling. 

Inspired by the above observations we now derive an 
analytical model for the jet formation as follows: As a 
flrst step the flow fleld of the collapsing cavity before 
pinch-off will be described by a line of sinks along the 
axis of symmetry Q . The strength of these sinks will 
be determined from the simulation at pinch- ojf and forms 
the only input quantity required by our model. Next, we 
will show how this picture naturally leads to a good de- 
scription of the bulk flow after pinch-off. The line of sinks 
acquires a "hole" between the two jets and an additional 
point sink emerges near the jet bases. Finally, we will 
obtain two differential equations for the widening and 
upward motion of the jet base which are the two most 
relevant processes for jet formation. Secondary processes 
as jet breakup and the precise dynamics of the jet tip are 
not addressed in this work. 

As a starting point. Green's identity allows us to write 
the potential at any point r in the liquid bulk as an in- 
tegral of sources and dipoles over the free surface: 



47r0 (r) = / d5'n' 
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with the integration taken over the free surface S as il- 
lustrated in Fig. [5] (a) and (b). Since the dipole term 
decays quickly as 1/ |r — r'|^, the source term (which de- 
cays only as 1/ |r — r'|) will be the only relevant contri- 
bution to the integral if the observation point is chosen 
sufficiently far from the free surface. As the cavity close 
to pinch-off becomes slender, d(j)/dn w for a point R 
on the free surface. Since the surface has no overhangs we 
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FIG. 4: (a) Sketch of the coUapsing cavity being described by 
a distribution of sinks (orange arrows) on its free surface, (b) 
During jet formation the cavity collapse in the outer region 
remains unchanged (orange arrows) while around the jet base 
sinks accumulate (green arrows). This can be approximated 
as a line of sinks along the axis of symmetry plus a point sink 
(green dot). In the central region around the pinch-off point, 
a hole is formed: sinks are completely absent. For a detailed 
description see main text. 



write dS = 2nRdz. Approximating the radial distance 
as r — r' « r turns Eq. (jT)) into [y] : 

2cP{r,z)= [ -^^^k^2^=dz' . (2) 
^ y^r2 + (z - z'f 

with a time- and height-dependent hne distribution of 
sinks qa.xis{z,t) along the axis of symmetry. Keeping in 
mind the extremely short time scale of jet formation as 
compared to the cavity collapse, we can assume the sink 
strength to remain constant in time from the moment of 
pinch-off tc onwards, qdz) = gaxis(2,ic)- 

During jet formation we divide the free surface into 
two regions separated by the jet base. The outer region 
contains the collapsing cavity until the jet base, while 
the inner region extends from the base inwards to the 
axis of symmetry as sketched in Fig. [4] (b). The princi- 
pal fluid motion in the outer region remains identical to 
the collapsing cavity before pinch-off. High up in the jet, 
the motion will be vertically upwards with negligible ra- 
dial velocity and thus will not contribute to the integral 
Eq. (21). As a free surface particle travels through the jet 
base and further up into the jet, it transitions from one 
flow regime to the other by decelerating its initial radial 
motion and turning it into vertical momentum. Thereby, 
its contribution to the integral ^ decays to a negligible 
amount. This decay of the sink strength cannot hap- 
pen instantaneously which leads to an accumulation (see 
Fig. S]) of sinks around the jet base and a correspond- 
ing inward motion in that area. The length over which 
the sinks decay and accumulate must be proportional to 
the radius of the jet base which is the only relevant local 
length scale, Cri,, with C a constant of order one. 

At this point it is good to stress that our model is 
not constructed to describe the liquid flow inside the jet 
itself but is valid only for the bulk flow outside the actual 



jet. This is sufficient for the explanation of jet formation 
which is the focus of the present paper. The sinks on 
the axis thus always remain outside of the liquid domain 
which they aim to describe. 

From an observation point located at r ^ ri,, the 
contribution of the sinks accumulating around the base 
can be regarded as a point sink of strength Crbqc{zb) 
since qdz) ~ qdzb) along the length Crf,. The point 
sink will be located a small distance above the base 
which is again proportional to the local length scale, i.e., 
Zs'mk = Zb + Csink ' fb introducing a second constant Csink 
of order unity. 

Similarly, the most relevant contribution of the outer 
region will be that part of the line integral closest to 
the observation point r. For an observation point at an 
altitude similar to or lower than the jet base, this is the 
region close to the jet base where again qdz) ~ qd^b)- 
To allow analytical treatment of the integral resulting 
from Eq. we can thus at any given time assume a 
sink strength being constant in space along the entire 
axis above the jet base. Through the motion of the jet 
base this sink strength depends implicitly on time qb{t) = 

qdzb{t)). 

Combining the approximations of the preceding para- 
graphs, we are now able to give an analytical expression 
derived from Eq. ([2]) for the potential at any point (r, z) 
as a function of the base position and z^: 

r°° dz' dz' 

' i-oo + {z- z'f ' L., + (z - z'f 

" V ' V ' 

collapsing cavity hole 



1/7-2 + (z - (Zfc -f- CsinkT'b))' 




point sink 

The initial sink distribution is obtained from the numer- 
ics by calculating qdz) = —RR along the surface just 
once at pinch-off. It forms the only input quantity re- 
quired by our jetting model. Note that, as we are dealing 
with the upwards jet, the effect of the point sink for the 
downward jet can be neglected. 

In order to derive the desired ODEs for ri,{t) and 
Zb(t) wc apply the Bernoulli equation with zero pressure 
d(j)/dt + \ V4>\ /2 = on the free surface. Wc then em- 
ploy Eq. ([3]) to obtain the first differential equation in- 
volving rb{t) and Zb{t). The second ODE results from 
the application of the kinematic boundary condition at 
the jet base d(t>/dz = dz^/dt. This leads to a closed sys- 
tem of two ODEs. The calculations are straightforward 
but lengthy and are omitted here. With C — 4.55 and 
C'sink = 0.63 the agreement with simulations and experi- 
ment is remarkable as demonstrated by Fig. [21 Wc stress 
that the model requires as its only ingredient the sink 
strength distribution at pinch-off. 

Finally, it is important to understand which region of 
the liquid bulk at pinch-off will eventually become ejected 
into the jet. This knowledge can be obtained from the 
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FIG. 5: A line of tracer particles is injected at the base of 
the jet at t/To = 0.08 (a). After advecting them backwards 
in time to the moment of pinch-off they form the border of a 
thin surface layer containing the liquid that eventually ends 
up in the jet (b). A movie in which tracers are inserted at the 
base of the jet and advected backward in time is included as 
EPAPS document No. XXXXXX. The thickness of this layer 
is plotted in (c). 

numerical simulations by injecting a line of particles at 
the base of the jet, cf. Fig. [5] (a). Since the flow field 
is known for all times previous to particle injection, the 
tracers can be followed backwards to their origin at i = 0. 
The line of tracers injected at the final instant will yield 
the outer boundary of the fiuid layer that, together with 
the free surface, delimitates the fiuid volume from which 
the jet originates. While the radial extent of the fluid 
layer depicted in Fig. [5] (c) is of the order of the disc ra- 
dius, its maximum thickness is only about 0.01 -iJo- Thus 
far, a similar surface layer has only been observed when 
jetting is directly caused by surface waves In the 

present case, the thinness of the layer is even more re- 
markable as it does not arise from a surface phenomenon 



but from the collapsing motion of the entire bulk liquid. 

In conclusion, we have studied in detail the mechanism 
responsible for the formation of high-speed Worthington 
jets after the impact of solid objects on a liquid surface. 
We showed that the liquid forming the jet originates from 
a thin layer straddling the surface of the impact cavity. 
Our main finding, nevertheless, is the vital importance of 
the radial energy focussing along the entire wall of this 
cavity. In contrast to other situations 0, H Q , the hyper- 
bolic flow around the singular pinch-off point turned out 
to be not the relevant mechanism behind jet formation. 
Instead, our case seems more reminiscent of the violent 
jets observed during the explosion of lined cavities [2l[. 
We proposed an analytical model which is in very good 
quantitative agreement with experimental data and nu- 
merical simulations. The only ingredients to the model 
are two constants of order one and a sink distribution 
qc{z) describing the collapsing cavity at pinch- off. 

We expect that the present mechanism is also respon- 
sible for the very thin jets ejected after the impact of 
water droplets on a liquid pool in a parameter range 
where a small cylindrical cavity at the bottom of the 
crater collapses in a very similar fashion as the impact 
cavity described in this work. In the future, our model 
of jet formation can serve as the base for predicting the 
shape and the velocity of the jet itself. 
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